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A parameter study is underway in Heidelberg that aims to find and constrain possible 
solutions to the dSph satellite problem under the assumption that some of these 
systems may not be dark matter dominated. The present state of the parameter 
survey of tidally disrupting spheroidal dwarf galaxies is described, and examples of 
preliminary results are presented. 
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1. Introduction 



At least about ten dwarf spheroidal (dSph) galaxies are known to orbit the Milky Way at 
J> ' distances ranging from a few tens to a few hundred kpc. On the sky they are barely discernible 

stellar density enhancements. Their velocity dispersions are similar to those seen in globular 
, clusters and they have approximately the same stellar mass. However, they are about two 

orders of magnitude more extended. For spherical systems in virial equilibrium with an isotropic 
velocity dispersion, the overall mass of the system can be determined from the observed velocity 
dispersion. Comparing this 'gravitational' mass to the luminosity of the system determines the 
mass-to-light ratio, M/L (in the following always given in solar units Mq/Lq). Values for M/L 
of about 10 or larger are usually taken to imply the presence of dark matter in a stellar system. 
For the dSph satellites, M/L values as large as a few hundred are inferred, implying that these 
systems may be completely dark matter dominated (for reviews see Mateo 1997, 1998). 

An alternative possibility relying on Newtonian physics, is that some of the observed dSph 
satellites may not be in virial equilibrium and may not be spherical but with non-isotropic 
velocity dispersions, in which case the observed (M/L) Q b s ratios may not be physical. That such 
systems may exist is shown by the simulations of the long-term evolution of initially spherical 
satellite galaxies with masses of 10 M@, initial Plummer radii of 300 pc on eccentric orbits in 
an extended Galactic dark halo with a circular velocity of 200 km/s (Kroupa 1997, hereinafter 
K97). The finding is that after about 99 per cent of the mass is lost from the satellites through 
repetitive tidal modification, remnants with (M/L) Q \, S > 10 remain that survive for a few orbital 
periods. This study is extended by Klessen & Kroupa (1998, hereinafter KK98) covering more 
orbits and Galactic halo masses using two different numerical codes. 

A shortcoming common to these remnants is that they have a central surface brightness 
too faint by about one order of magnitude, if (M/L) tme = 3 is the mass-to-light ratio of the 
stellar population. One possible solution to this problem is to postulate that (M/L) tme w 0.3 
(K97). However, this implies an unusual initial mass function. Another possibility is to scale 
the satellite up in mass and size so that each particle carries a larger mass (for simulations with 
the same number of particles) and thus light, while keeping the binding energy of the satellite 
roughly unaffected (KK98). This, however, leads to remnants that have too large half-light radii. 

It is thus clear that the simulations mentioned above, as well as other observational evidence 
summarised in K97 and KK98, give a strong hint that some of the Galactic dSph satellites may 
not be dark matter dominated but that they may be significantly affected by tides. However, 
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much work remains to be done to establish if models can be produced that fit all the observational 
properties of at least some of the known dSph satellites. 

It is the purpose of this contribution to give a preliminary account of the extensive parameter 
survey under way in Heidelberg that aims at constraining the possible region in (initial satellite 
binding-energy, mass and extend of the Galactic dark matter halo, and orbit) parameter space 
in which viable "dSph solutions" exist in the tidal scenario. If none such region can be found 
then the conclusion that all of the dSph satellites are dark matter dominated is strengthened. 

2. Simulations 

The study of parameter space is performed with Superbox, which is a particle mesh code with 
nested sub-grids. It is described by Madejski & Bien (1993) with a brief description of the 
modern version given by K97. A more thorough account will be available in Fellhauer et al. 
(1998). The reason for using Superbox lies in its tremendous efficiency: the self-consistent 
interaction of two galaxies consisting in total of 1 — 2 x 10 6 particles can be computed on small 
workstations and personal computers with Pentium processors within a few days. 

The problem under investigation here is reduced to the interaction of a low-mass satellite 
galaxy with a massive Galactic dark halo. The details of setting up the satellite and Galactic halo 
models can be found in K97. The dark halo is assumed to be an isothermal mass distribution 
with a circular velocity of 200 km/s and a cutoff radius of R c . The satellite is defined by its 
initial mass, M sat , and its initial Plummer radius, Rp\. The satellite is placed into the Galactic 
dark halo at apo-galactic distance -R apo and with a tangential velocity t> t . The resulting orbit 
has an eccentricity e = (Ra, po — Rper\) / {R&po + Rperi), where i? per i is the peri-galactic distance. 

The satellite is observed as an observer from Earth sees a dSph satellite. It's projected 
brightness profile is measured. The central surface brightness and half-light radius are computed 
through profile fitting, and the line-of sight velocity dispersion is determined. The hypothetical 
observer thus estimates (M/L) Q b s as a function of time. 

3. Some Results 

The parameter study proceeds by varying M sat , Rpi, R c , i? apo and v t - Most of the simulations 
done so far concentrate on a satellite with M sat = 10 7 M Q and R Pi = 300 pc in a Galactic dark 
halo with R c = 250 kpc or 40 kpc. 

Typically, (M/L) Q b s remains constant at its initial low value until disruption time, when the 
satellite looses most of its remaining mass during passage through peri-galacticon. Thereafter, 
a remnant with roughly 1 per cent of the initial satellite mass remains. It has properties rather 
similar to typical dSph satellites, with large (M/L) D b s . Evolution curves for (M/L) Q b s (i) and 
other satellite parameters are presented in K97 and KK98. 

A particularly noteworthy result obtained by KK98 is that a remnant has a line-of-sight 
velocity dispersion a > 6 km/s if e > 0.5. This suggests a useful observational criterion, for if 
e < 0.5 for a > 6 km/s then the respective system is probably dark matter dominated. Here, 
preliminary additional results are presented for this satellite, with some discussion of other 
models as well. 

The time when (M/L) Q b s > 50 occurs is denoted by T50; it increases for increasing e, as is 
evident from Fig. [I]. From the figure it is apparent that T50 is primarily a function of e and R apo - 
Two observational constraints are plotted: The dSph satellite Ursa Minor contains only old 
stars with ages r ~ 15 Gyr (van den Bergh 1994). It's preliminary proper motion (without an 
uncertainty) is reported by Olszewski (1998), from which the orbital eccentricity is estimated. 
The dSph satellite Sculptor has stars with ages in the range r ~ 5 — 15 Gyr (van den Bergh 
1994), and the proper motion has been estimated by Schweitzer et al. (1995). Both satellites 
are at a distance of 60-70 kpc. It is interesting to note that UMi has (M/L) ^ s = 95 ± 43, while 
Scl has (M/L) ohs = 10.9 ± 7.5 (Irwin & Hatzidimitriou 1995). 
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Figure 1: The time t = T50, when (M/L) Q b s (t) > 50, is plotted against orbital eccentricity for a 
number of simulations with M sa t = 10 7 Mq and Rp\ = 300 pc, as defined in the key. 



The application of Fig. |l] is as follows: If the entire satellite was assembled a few Gyr ago 
then this would imply some recent merger event strong enough to lead to the formation of tidal 
tail dwarfs (e.g. Kroupa 1998). It is not clear if there is any other evidence for such an event. 
The present masses of the stellar components (assuming (M/L)t rue = 3) of UMi and Scl are, 
respectively, about 6 x 10 5 Mq and 4 x 10 6 Mq (Irwin & Hatzidimitriou 1995), which may be too 
low to accrete gas from a hypothetical gas cloud on a similar orbit (c.f. Magellanic or similar 
stream), or let alone retain gas over time-spans of many Gyr. Any stars with nuclear ages of a 
few Gyr were probably born from some gas-accretion event during a more massive past. Perhaps 
one may identify the nuclear age, r, of the youngest stars as being the time when the satellite 
was massive enough for the last time to accrete some gas. The precursor satellite galaxies will 
have been two orders of magnitudes more massive than the dSph satellites, if the tidal scenario is 
applicable. Thus, T50 < r if the respective dSph satellite is in the remnant phase with artificially 
inflated (M/L) ohs . 

Unfortunately the observational constraints are too weak to allow any definitive conclusions. 
But it is clear that if a dSph with (M/L) Q b s > 10 containing relatively young stars (a few Gyr 
old, e.g. Scl)) can be shown with high confidence to be on a low eccentricity (e < 0.3) orbit, 
then the case for dark matter in this system would be rather strong. 

Ultimately, the parameter study will allow constraining that region in (M sa t,i?pi) space 
that may lead, within a Hubble time, to a remnant phase that agrees with the observational 
properties of at least some of the dSph satellites. Fig. ^ demonstrates this schematically. In this 
figure, the thick (D = 100 kpc) and thin (D = 50 kpc) short-dashed lines are the tidal radii, 
logio^t = — 3.51 + (l/3)log 10 M sa t + (2/3)log 10 -D in kpc, where M sat is in Mq and D is the distance 
from the Galactic centre in kpc. Satellites that lie initially to the right of these (Rp\ > r t ) are 
more or less immediately destroyed on injection into the dark halo. The long-dashed line, 
logioM sat = 8.3 + (l/2)log 10 -Rpi (units as above), is the binding energy, = (GM s 2 at )/i? p i, of 
the satellite. Satellites that initially lie above this line do not enter the remnant phase within 
a Hubble time. The filled dot is the satellite mentioned above, and the open circle represents a 
series of simulations that barely lead to any solutions, irrespective of e. 

Once the solution region is delineated in this diagram, inferences can be drawn about the 
physical properties of the progenitors of those dSph satellites that are not dark matter dominated. 
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Figure 2: Initial satellite mass as a function of initial characteristic satellite radius. 



This in turn may allow conclusions about the physical conditions that lead to the birth of such 
systems. 



4. Conclusions 

The simulations towards the parameter study briefly reported here have been running for a 
couple of years, having consumed about 1 CPU yr, and a total of more than 13 GB of data have 
accumulated. The present results do not unambiguously demonstrate that dSph-like systems 
without dark matter exist. Nevertheless, some of the models bear surprising similarities with 
the dSph satellites, and some inferences can already be drawn from this survey. 

There exists in (M sa t,i?pi) space a region that leads within a Hubble time to remnants of 
satellites that appear tantalisingly similar to the dSph satellites (K97). Predictions of correla- 
tions between observable parameters can be arrived at (KK98, Fig. |l|). Also, information as to 
the possible initial state of the progenitors of those dSph satellites that may not be dark matter 
dominated is emerging (Fig. ^). 
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